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Давление сыпучего материала в круглом силосе 
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Рассмотрена задача о распределении давления в сыпучем материале, находящемся в круглом силосе, с по- 
зиций вязкой жидкости, обладающей сухим внутренним трением. 
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Введение. Как известно, сыпучие материалы обладают как свойствами твёрдого тела, так и 
свойствами жидкости, что сильно осложняет теоретические исследования. Например, достаточно 
элементарная задача расчёта давления в сыпучем материале, находящемся в цилиндрической 
ёмкости, с позиций сплошного твёрдого тела [1] приводит к резкому расхождению с эксперимен- 
тальными данными. 

Та же задача, рассматриваемая с позиций дискретного твёрдого тела, решается в первом 
приближении только для частиц сферической формы, «уложенных» горизонтальными слоями с 
правильной геометрической «упаковкой» [2]. При этом учитываются коэффициент трения между 
сферическими частицами (коэффициент внутреннего трения |;) и коэффициент трения частиц о 
стены ёмкости (коэффициент внешнего трения |), которые определяются экспериментально. И 
хотя согласие с экспериментом здесь гораздо лучше, такой подход также вызывает справедливые 
возражения. 

Предпринята попытка решить сформулированную выше задачу, представляя сыпучий ма- 
териал в виде вязкой жидкости с теми же коэффициентами внутреннего (р) и внешнего (ре) тре- 
ния. При таком подходе нет необходимости задаваться формой частиц, их упругостью, упаковкой, 
размерами ит. д. 

Теоретический расчёт давления в круглом силосе. Рассмотрим сыпучий материал в глубо- 
ком круглом силосе (см. рис. 1). 
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Рис. 1. Сыпучий материал в круглом силосе и цилиндрическая система координат: Н — высота силоса; К — радиус силоса; 
2, р — текущие координаты элементарного кольца 


В силу симметрии задачи воспользуемся цилиндрической системой координат, в которой 
можно пренебречь зависимостью от азимутального угла, как показано на рис. 1. В слое сыпучего 
материала выделим элементарное кольцо, объём которого будет равным: 2.п.р. ф. ай. 

В стационарном случае сумма сил, действующих на элементарное кольцо, равна нулю. В 
проекции на ось О7 это условие имеет вид: 
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0=у:а.2.-п.р. ф.@-2.п:р.:@ф.а2: 





(1) 


где у — плотность сыпучего материала; Р, — вертикальное, параллельное силе тяжести давление 
в сыпучем материале; Р, — горизонтальное, перпендикулярное силе тяжести давление в сы- 
пучем материале. 

Здесь мы, в согласии с экспериментом, полагаем, что закон Паскаля для сыпучих тел не 
справедлив. 

В уравнении (1) первое слагаемое — сила тяжести элементарного кольца, второе — сила 
сопротивления движению за счёт градиента вертикального давления, третье — сила внутреннего 
сухого трения. 

В проекции на горизонтальное направление это условие имеет вид: 





(2) 


Здесь первое слагаемое — сила сопротивления сжатию (или растяжению) элементарного 
кольца, обусловленная градиентом горизонтального давления, второе — сила трения на горизон- 
тальных поверхностях кольца. 

Из уравнения (1), учитывая (2), получим дифференциальное уравнение для определения 
вертикальной компоненты давления: 

ОР, 
0=у-9.2-п.р.ф- 42 -2.п.р.ф.@. (1+4). —^ т (3) 
Проинтегрируем уравнение (3) по координате р в от 0 до ДА: 
0=п-у.9.Е?. 4 -п. 2. . (+0). 52 оС. (4) 


Поскольку при таком интегрировании элементарный о переходит в элементарный 
слой, то в последнем уравнении первое слагаемое представляет собой силу тяжести, действую- 
щую на единицу длины (вдоль радиальной оси) элементарного слоя толщиной с, второе слагае- 
мое — результирующую силу сопротивления (за счёт градиента давления и силы внутреннего 
трения), действующую на этот элементарный слой. Поскольку уравнение (4) справедливо и на 
поверхности р = А, то очевидно, что константа С представляет собой силу трения, действующую 
на боковую поверхность элементарного слоя со стороны стенки силоса и направленную против 
силы тяжести: 

СЕ: В 2 П.В чо ==р: Ко, +2 ПКЕ. (5) 


Здесь использована связь между вертикальной и горизонтальной компонентами давления 
через боковой коэффициент К: Р, =К.Р, (такая зависимость установлена в [3] и рекомендована 


в СНиП 2.02.02-85). 
Общий интеграл уравнения (4) с учётом (5) можно представить в виде: 


_у:9.К +С -еф +2 


а те 

2.ы..К (1+7). ^ 

Если найти константу С, из условия, что на поверхности силоса давление равно нулю, то 

выражение для вертикального давления (без учёта зависимости от горизонтальной координаты) 
примет вид: 


(6) 





В = у 9 К 1-ехр = 


ие (+12). ий 
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Соответственно, выражение для горизонтальной компоненты давления запишется в виде: 


а | ехр 2.н..К.2 
2.Н. (1+52).^ 





(8) 


Для анализа полученного решения найдём разность Д между весом сыпучего материала в 
силосе и суммой сил давления на дно силоса Ар и силы трения сыпучего материала о боковую 
стену Ар. Нетрудно убедиться, что эта разность будет равна: 


_п.у. 9.3 5 , 2..-К.Н 
р 2.[.-К (1+ы и льет (1+5). ^ 








(9) 

Из данного выражения видно, что ЛД =0 в том случае, когда р; = 0. Таким образом, при 
выполнении условия равновесия движение сыпучего материала не происходит и внутренним тре- 
нием можно пренебречь. Из (9) также следует, что если р; + 0, то сила трения в сумме с силой 
давления будут меньше силы тяжести, за счёт чего, собственно, и происходит уплотнение сыпу- 
чего материала в силосе. 

Очевидно, что при р; -› 0, К-» 1 наступит равновесное состояние (сыпучее тело перехо- 
дит в идеальную жидкость), давление по горизонтали станет равным давлению по вертикали и, 
соответственно, решение (7) примет вид: 


у. 9.К 2-4. .2 
М ЖЕ. о 


е 





Графики зависимости давлений, рассчитанных по формулам (7) и (10), приведены на рис. 2. 
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Рис. 2. Вертикальное (сплошная линия) и горизонтальное (пунктирная линия) давление в стационарном состоянии и 
давление (линия из точек) в равновесном состоянии для круглого силоса с сыпучим материалом (К = 3,5 м, Х = 300°, 
у= 800 кг/м?) 


Из рисунка следует, что вертикальное давление в стационарном состоянии (сплошная ли- 
ния) больше давления в равновесном состоянии (линия из точек). В то же время горизонтальное 
давление в стационарном состоянии меньше давления в равновесном состоянии (пунктирная ли- 
ния). Следовательно, в равновесном состоянии стены силоса берут на себя большую нагрузку, в 
результате чего давление на дно уменьшается. Иначе говоря, в стационарном состоянии (при 
движении сыпучего материала) на стены силоса оказывается меньшая нагрузка, чем в равновес- 
ном состоянии, когда сыпучий материал неподвижен. 

При нахождении сыпучего материала в силосе формально он не перемещается, следова- 
тельно, находится в равновесном состоянии. При этом давление в нём необходимо рассчитывать 
по формуле (10). Но в действительности в сыпучем материале скачкообразно происходят процес- 
сы его уплотнения, слёживания под действием микроколебаний сейсмического или техногенного 
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характера, т. е. он находится в состоянии микродвижения. Таким образом, сыпучий материал 
часть времени «проводит» в стационарном состоянии, а часть — в равновесном. Соответственно, 
давление также меняется скачкообразно. 

Введём теперь в решение (7) зависимость от координаты р. Аналогично [3] эту зависи- 
мость представим в виде: 


узо 2...К 2-(К-р).19Х 
В = .|1 В , 
2 К ы 1+ [3 (11) 





е 
где х — угол естественного откоса. 
Потребуем теперь, чтобы решение (11) удовлетворяло уравнению (2) с учётом того, что 
Р.=К.Р: 
р 2 








ко 98 ь 
ОЕ о? (12) 


Найдя первые производные от (11) по координатам р и ги подставив их в равенство (12), 
получаем, что оно выполняется только при К. 9х =Ы, ‚т. е. 


Ы, 
ы. 13 
9х (13) 


Это выражение для бокового коэффициента совпадает с результатом, полученным в [3] 
для сыпучего материала, находящегося в насыпи. 

Однако решение в форме (11) не удовлетворяет граничному условию на поверхности 
= 0. Поэтому для нахождения более точного решения разобьём весь объём силоса, занятый сы- 
пучим материалом, на области (рис. 3). 











Рис. 3. Деление силоса на области: Н — высота силоса; А — радиус силоса; х — угол естественного откоса; 2, р — текущие 
координаты элементарного кольца; Аь — произвольная точка на границе областей Ти П 


Поскольку в области 1 влияние стенок силоса ещё не будет сказываться (эта область рас- 
положена выше плоскости естественного откоса), то в этой области решение будет определяться 


формулой [3]: 
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Р!(2)= "9.2. 
а И. (14) 
Решение в области 1 будем искать в виде: 
ри ря зу 9 .|1-ех е . + . (Ю .Ё ; 
‚ (р,2) та а г. ти р) 9х (15) 


На линии 2 = 0 решение (15) обращается в нуль, на боковой стенке силоса (при р = А) оно 
переходит в уравнение (7), а на поверхности естественного откоса (на границе областей), опре- 
деляемой уравнениями р =А,, 2 = (к -К,)- 9х ‚ уравнение (15) переходит в уравнение (14). Та- 
ким образом, согласование решений в областях Г и П выполнено. 

В случае, когда сыпучий материал находится в силосе с абсолютно гладкими стенами 
(н. >00), решение (15) примет вид: 








ит Ри = т :|[2- (В р). 9х] +-^ Я. (к -р)- 9х = ть и. 
> 1+Н; 1+. 1+ы; 

Это означает, что сыпучий материал в этом случае ведёт себя как безграничная насыпь и 
влиянием стенок можно пренебречь, что не противоречит физическому смыслу задачи. 

Итак, мы видим, что решения (14) и (15) удовлетворяют граничным условиям и предель- 
ному переходу к идеальной жидкости. 

Таким образом, для идеального сыпучего материала в круглом силосе найдено вертикаль- 
ное давление (15). Горизонтальное давление определяем из выражения (15) с учётом выражения 
(13) для бокового коэффициента. 

На рис. 4 показана зависимость вертикального и горизонтального давления, рассчитанно- 
го по формулам (14) и (15), в силосе радиуса А = 3,5 м, в котором находится сыпучий материал 
(\ = 800 кг/м3, Х = 30°), непосредственно вдоль его стенки и по центру силоса. Как видно из ри- 
сунка, вертикальное давление по осевой линии силоса (р = 0) несколько выше вертикального 
давления у стенки силоса (р = А) и, кроме того, горизонтальное давление меньше вертикального. 





область 1 


область П 





0 5 10 15 р. > 
2 м 


Рис. 4. Зависимость вертикального и горизонтального давления от координаты 2: Р, (0, 2) — вертикальное давление на оси 
силоса; Р, (К, 2) — вертикальное давление у стенки силоса; Р, (К, 2) — горизонтальное давление у стенки силоса; радиус 
силоса А = 3,5 м; угол естественного откоса х = 30°; плотность сыпучего материала у = 800 кг/м?; р; = 0,287; ие = 0,364 


На рис. 5, 6 представлены графики сравнения теоретических кривых (непрерывные ли- 
нии) и экспериментальных значений [4]. На рис. 5 показана зависимость горизонтального давле- 
ния на стену круглого железобетонного силоса высотой Н = 21 м и радиусом К = 3,25 м, а на 
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рис. 6 — зависимость силы давления на дно стального круглого силоса высотой Н = 6 ми радиу- 
сом А = 0,5м. 
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Рис. 5. Теоретическое горизонтальное стационарное (сплошная линия) и равновесное (линия из точек) давление в 
круглом железобетонном силосе (Н = 21 м, К = 3,25 м, Хх = 30°); вертикальные линии — разброс экспериментальных 
значений 
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Рис. 6. Теоретическая сила давления на дно стального круглого силоса (Н = 6м, А = 0,5 м) от высоты засыпки: сплошная 
линия — стационарное состояние, линия из точек — равновесное (х = 30°); вертикальные линии — разброс 
экспериментальных значений 
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Рис. 7. Зависимость силы трения на боковую стенку круглого стального силоса (Н = 6 м, К = 0,5 м, Х = 30°) от высоты 
засыпки: стационарное состояние — сплошная линия, равновесное — линия из точек, вертикальные линии — разброс 
экспериментальных значений 
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На рис. 7 представлена зависимость силы трения, действующей на боковую стенку сталь- 
ного круглого силоса высотой Н = 6 ми радиусом К = 0,5 м, от высоты засыпки. 

Как видно из рисунков, разброс экспериментальных данных достаточно велик. Кроме того, 
само горизонтальное давление по направлениям сторон света разное [4]. 
Выводы. Очевидно, что этот разброс экспериментальных значений связан с процессами уплот- 
нения сыпучего материала, со временем его выдержки в силосе, а также с тем, с какого направ- 
ления идут микровибрации сейсмического или техногенного характера. Сыпучий материал в сило- 
се всё время находится в промежуточном положении между стационарным и равновесным со- 
стояниями. Кроме того, сами параметры сыпучего материала (насыпная плотность, углы внешнего 
и внутреннего трения) являются некоторыми обобщёнными, мало соответствующими конкретному 
сыпучему материалу и, к тому же, меняющимися с течением времени. Поэтому более точного со- 
гласия с экспериментальными данными ожидать трудно. 

Но из сравнения с экспериментальными данными всё же можно сделать вывод, что со- 
стояние сыпучего материала ближе к стационарному, чем к равновесному. 

Таким образом, подход к сыпучему материалу как к вязкой жидкости приводит к согла- 
сию с экспериментальными данными при минимальном количестве необходимых допущений. 
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